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^s^j , Anti-decuplet penta-quark baryon is studied with the quenched anisotropic lattice 

f— *^ ■ QCD for accurate measurement of the correlator. Both the positive and negative 

f*^ ' parity states are studied using a non-NK type interpolating field with 1 = and 

J = 1/2. After the chiral extrapolation, the lowest positive parity state is found 
at rrtg, ~ 2.25 GeV, which is too massive to be identified with the experimentally 
i observed + (154O). The lowest negative parity state is found at me ~ 1.75 GeV, 

"■»«. ' which is rather close to the empirical value. To confirm that this state is a compact 

5Q resonance, a new method with "hybrid boundary condition (HBC)" is proposed. 
The HBC analysis shows that the observed state in the negative parity channel is 
JL ■ an NK scattering state. 

M-t, 

1. Introduction 

LEPS group at SPring-8 has discovered a narrow resonance O + (1540), 
which is centered at 1.54 ± 0.01 GeV with a width smaller than 25 MeV. 1 
This resonance is confirmed to have baryon number B = 1, charge Q = +1 
and strangeness S = +1 implying that it is a baryon containing at least 
one s. Hence, its simplest configuration is uudds, i.e., a manifestly exotic 
penta-quark (5Q) state. The experimental discovery of + was motivated 
by a theoretical prediction. 2 

Tremendous theoretical efforts have been and are still being devoted 
to the investigation of 8 + , 3 ' 4 among which its parity is one of the most 
important topics. Experimentally, the parity determination of + is quite 
challenging, 5,6 while opinions are divided in the theoretical side. 3 

There are several quenched lattice QCD studies of the 5Q state, 7 ' 8,9 ' 10 
However, the results have not yet reached a consensus. One group 9 claims 
the existence of a low-lying positive parity 5Q resonance. Negative parity 
5Q resonance is claimed by two groups, 7,8 among which Ref. 8 has omit- 
ted a quark-exchange diagram between diquark pairs assuming the highly 
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'Lattice QCD numerical calculation has been done with NEC SX-5 at Osaka University. 
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correlated diquark picture. Note that these three groups employed non-NK 
type interpolating fields. In contrast, Ref. 10 has employed the NK-type 
interpolating field, and performed solid analysis concluding that no signal 
for a 5Q resonance state is observed. There is another type of lattice QCD 
studies of the static 5Q potential 11 aiming at providing physical insights 
into the structure of 5Q baryons. 

In this paper, we study the 5Q baryon 6 + for both parities by using 
high-precision data generated with the quenched anisotropic lattice QCD. 
We employ the standard Wilson gauge action at [3 = 5.75 on the 12 3 x 96 
lattice with the renormalized anisotropy a s /a t = 4. The anisotropic lattice 
method is a powerful technique, which can provide us with high-precision 
data quite efficiently. 12 ' 13,14,15 The lattice spacing is determined from the 
static quark potential adopting the Sommer parameter r^ 1 = 395 MeV 
leading to a" 1 = 1.100(6) GeV (a s ~ 0.18 fm). 13 The lattice size 12 3 x 96 
amounts to (2.15fm) 3 x 4.30fm in the physical unit. The (9(a)-improved 
Wilson quark (clover) action is employed 13 with four values of hopping 
parameters as k = 0.1210(0.0010)0.1240, which roughly covers m s < m q < 
2m s corresponding to m v /m p = 0.81, 0.77, 0.72 and 0.65. By keeping k s — 
0.1240 fixed for s quark, we change n = 0.1210 — 0.1240 for u and d quarks 
for chiral extrapolation. Unless otherwise indicated, we use 



as a typical set of hopping parameters. Anti-periodic boundary condition 
(BC) is imposed on the quark fields along the temporal direction. To en- 
hance the low-lying spectra, we adopt a smeared source with the gaussian 
size p ~ 0.4 fm. We use 504 gauge configurations to construct correlators 
of + . For detail, see Ref. 16. 

In the former part of this paper, we present the standard analysis of 5Q 
correlators in both the positive and the negative parity channels adopting 
the standard periodic boundary condition along spatial directions. Latter 
half of this paper is devoted to a further investigation of the negative parity 
state. Proposing a new general method with "hybrid boundary condition 
(HBC)" , we attempt to determine whether it is a compact resonance state 
or a NK scattering state. 

2. Parity projection 

We consider a non-NK type interpolating field for 9+ as 



{k s ,k) = (0.1240,0.1220), 



(1) 



O = e abc e ade e bfg (u^C-f 5 d e ) (ujCdg) (CsJ) 



(2) 
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where a — g denote color indices, and C = 7472 denotes the charge conjuga- 
tion matrix. The quantum number of O is spin J = 1/2 and isospin 1 = 0. 
Under the spatial reflection of the quark fields, i.e., q(t,x) — ► j4q(t 7 — x), 
O transforms exactly in the same way, i.e., 0(t,x) — > +"f^O(t, — x), which 
means that the intrinsic parity of O is positive. Although its intrinsic parity 
is positive, it couples to negative parity states as well. 17 

We consider the asymptotic behavior of the correlator in the 5Q CM 
frame as 

G a/3 (t) = i^O a (t, 5)0^(0,5)), (3) 

X 

where V denotes the spatial volume. In the region of -C t -C N t (N t : 
temporal lattice size), the correlator is decomposed into two parts as 

G(t) = P+ (c+e- m+t - C-e-" 1 -^*-*^ (4) 

+ P- {C-e- m - 1 - C+e-"^^-^ , 

where m± refer to the energies of lowest- lying states in positive and negative 
parity channels, respectively P± = (1 ± ji)/2 serve as projection matrices 
onto the "upper" and "lower" Dirac subspaces, respectively, in the standard 
Dirac representation. Eq. (4) suggests that, in the region of -C t <C N t /2, 
the backwardly propagating states can be neglected. Hence, "upper" Dirac 
subspace is dominated by the lowest-lying positive parity state, whereas 
"lower" Dirac subspace is dominated by the lowest-lying negative parity 
state. We utilize this property in parity projection. 

3. Numerical Result with standard BC 

In Fig. 1, we show the effective mass plots for both the parity channels 
adopting Eq. (1). The effective mass is defined as 

m cB (t) = \og{g(t)/g(t + 1)), (5) 

where g(t) denotes correlator in Dirac "upper" or "lower" subspaces. For- 
mally, g(t) can be expressed as a sum of exponentials. In the asymptotic 
region <C t <C N t /2 in Euclidean time, contaminations of excited states 
are expected to be reduced. If g(t) is dominated by single exponential 
corresponding to the lowest-lying mass m, then the effective mass behaves 
as constant in this region, i.e., m e s(t) ~ m. Owing to this property, the 
effective mass plot is often used to determine the fit range. 17 
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For both the parity channels, we find plateaus in the region 25 < t < 
35, where single-exponential dominance is expected to be achieved. We 
simply neglect the data for t > 35, where contributions from the backward 
propagations are seen to become less negligible. The single-exponential 
fit is performed in the plateau region. The results are denoted by solid 
lines. The dotted lines indicate the p-wave (s-wave) NK thresholds for 
positive (negative) parity channels on the spatial lattice size L ~ 2.15 
fm. Note that due to the quantized spatial momentum in the finite box, 
the p-wave threshold is raised as E t h — \/ m % 
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Figure 1. The effective mass plots of positive and negative parity © + adopting Eq. (1). 
The solid lines denote the result of the single-exponential fit performed in the region, 
25 < t < 35. The dotted lines denote the p-wave (s-wave) NK threshold energy for 
positive (negative) parity channels on the spatial lattice size L ~ 2.15 fm. 



In Fig. 2, the masses of positive (triangle) and negative (circle) parity 
+ arc plotted against mj. The open symbols denote direct lattice data. 
We find that the data behaves linearly in mj. Such a linear behavior 
against is also observed for ordinary non-PS mesons and baryons. 13 ' 14 
We extrapolate the lattice data linearly to the physical quark mass region. 
The results are denoted by closed symbols. For convenience, we show p- 
wave (upper) and s-wave (lower) NK threshold with dotted lines. 

In the positive parity channel, the chiral extrapolation leads to me = 
2.25 GeV. Since it is too massive, it cannot be identified with the experi- 
mentally observed 6 + (1540). In contrast, in the negative parity channel, 
the chiral extrapolation leads to me = 1.75 GeV, which is rather close to 
the empirical value. In order for this to be identified with 9 + (1540), it 
should be confirmed that the observed state is not a NK scattering state 
but a compact 5Q resonance state. We will pursue this direction in the 
next section. 
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Figure 2. m@ for both parity channels against mj.. The triangles correspond to the 
positive parity, while the circles correspond to the negative parity. The open symbols 
denote direct lattice data, whereas the closed ones the results after the chiral extrapola- 
tion. The dotted lines indicate the NK threshold energies for p-wave (upper) and s-wavc 
(lower) cases. 

4. Further investigation with hybrid BC 

In the positive parity channel, NK scattering states are in p-wave. Hence, 
in the 5Q CM frame, the minimum momenta of N and K are non-zero, 
i.e., |p m ; n | = 2n/L (L : the spatial size of the lattice), through which the 
NK threshold energy acquires the explicit volume dependence as E t h — 
\J m 2 N + p^ in + \Jrr? K + pfjjin- This can be utilized to determine whether 
the state of concern is a compact resonance or a NK scattering state. On 
the other hand, in the negative parity channel, NK scattering states are 
in s-wave. Hence, in the 5Q CM frame, N and K can have zero spatial 
momentum, and the NK threshold energy does not have an explicit volume 
dependence, i.e., E t h ~ + tuk- This is not convenient for our purpose. 

We may ask ourselves whether there could be some prescription which 
can provide s-wave NK threshold energy with an explicit volume depen- 
dence as the p-wave one. This is achieved by twisting the spatial boundary 



Table 1. The hybrid boundary condition (HBC) imposed 
on the quark fields. The second line shows the standard 
(periodic) BC for comparison. 





u quark 


d quark 


s quark 


HBC 


anti-periodic 


anti-periodic 


periodic 


standard BC 


periodic 


periodic 


periodic 



condition (BC) of quark fields in a flavor dependent manner as follows. We 
impose the anti-periodic BC on u and d quark fields, while periodic BC on 
s quark field. We will refer to this boundary condition as "hybrid boundary 
condition (HBC)". (See Table 1.) 

Under HBC, hadrons are subject to their own spatial BC. For instance, 
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since N(uud,udd) and K(us, ds) contain odd numbers of u and d quarks, 
they are subject to the anti-periodic BC. In contrast, since <d + (uudds) 
contains even numbers of u and d quarks, it is subject to the periodic BC. 
(See Table 2.) Recall that, in the box of the size L, the spatial momenta 



Table 2. The consequence of the HBC on hadrons 





quark content 


spatial BC 


minimum momentum 




N 


uud,udd 


anti-periodic 


(±n/L,±ir/L,±n/L) 


IPminl = VStt/L 


K 


us,ds 


anti-periodic 


(±tt/L,±tv/L,±7t/L) 


|Pmin| = V37T/L 


e+ 


uudds 


periodic 


(0,0,0) 


IPminl = 



are quantized as pi = 2niTr/L for periodic BC and (2rij + l)7r/L for anti- 
periodic BC with rii e 2. Therefore, 9 + can have zero spatial momentum 
as |pmin| = 0. In contrast, N and K have non-zero minimum spatial mo- 
menta as |pmin| = VStt/L. Thus, under HBC, s-wave NK threshold energy 
is raised as 

E th - \Jm 2 N +p 2 min + ^Jm^+p 2 ^, |p min | = V3n/L, (6) 

whereas a compact 5Q resonance state is expected to be unaffected. 

In order to see that HBC does not affect the spatially localized resonance 
states, we show an example of an established resonance S(wds) baryon. We 
select S, because it is subject to the periodic BC unlike N. In Fig. 3, we show 

Standard BC Hybrid BC 

E (uds;j"=l/2~) 



' ■ 
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Figure 3. The effective mass plots of T,(uds) under the standard BC and HBC adopting 
Eq. (1). Solid lines denote the best-fit results performed in the region 20 < t < 30. 

the effective mass plots of S for the standard BC and HBC. The solid lines 
denote the best-fit results performed in the plateau region as 20 < t < 30. 
There is no significant difference between the two best-fit masses, which 
shows that HBC does not affect the localized resonance states. 

Now we present the HBC result of the 5Q effective mass in the negative 
parity channel in Fig. 4. The dotted line denotes the modified NK threshold 
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Figure 4. The effective mass plot for the negative parity © + under HBC adopting 
Eq. (1). The solid line denotes the result of the best-fit performed in the plateau region 
as 25 < t < 35. The dotted line denotes the s-wave NK threshold Eq. (6). This figure 
should be compared with the r.h.s. in Fig. 1. 

energy due to HBC. Note that the shift of the NK threshold amounts to 
about 200 MeV in this case, i.e., L ~ 2.15 fm with Eq. (1). We observe 
that the plateau is raised by consistent amount as the shift of the threshold, 
which shows that there is no compact 5Q resonance in the region as 



m N + m K <E < ^Jm 2 N + p^ in + \J m 2 K + p^ in , |p min | = V3n/L. (7) 

In particular, the plateau observed in the negative parity channel in the 
previous section turns out to be an NK scattering state. 

In Fig. 5, we show the comparison of the results of standard BC (l.h.s.) 
and HBC (r.h.s.) for each k. Dots denote the best-fit mass obtained in 
their plateau region, the solid lines denote the NK threshold energies. We 
see that, for all k, the plateaus are raised in a consistent amount as the 
shift of the NK threshold energies. 
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Figure 5. Comparison of the results of standard BC (l.h.s) and HBC (r.h.s.) for each k. 
Dots denote the best-fit mass obtained in their plateau regions. The solid lines denote 
the corresponding NK threshold energies. 
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5. Summary and Discussion 

We have studied the penta-quark 9 + state with the quenched anisotropic 
lattice QCD to provide high-precision data. After the chiral extrapolation, 
we have obtained m@ — 2.25 GeV in the positive parity channel. Since it is 
too massive, we have concluded that it cannot be identified with the exper- 
imentally observed 9 + (1540). In contrast, in the negative parity channel, 
we have obtained me = 1-75 GeV, which is rather close to the empiri- 
cal value. In order to confirm that this state is a compact 5Q resonance, 
we have proposed a new general method with "hybrid boundary condition 
(HBC)" . The HBC analysis has showed that the observed states in the 
negative parity channel are NK scattering states for all values of k. 

We have thus observed no relevant signals on the compact 5Q resonance 
in both the parity channels. To reveal the mysterious nature of 9 + (1540), 
more systematic investigations seem to be necessary. In the future study, it 
is desirable to examine the large volume effect, the dynamical quark effect, 
different interpolating fields including highly non-local ones, and different 
quantum numbers other than J = 1/2, 1 = 0. 

References 

1. LEPS Collaboration, T. Nakano et al, Phys. Rev. Lett. 91, 012002 (2003). 

2. D. Diakonov, V. Petrov and M.V. Polyakov, Z. Phys. A359, 305 (1997). 

3. For a review, M. Oka, Prog. Theor. Phys. Ill, 1 (2004) and its references. 

4. S.L. Zhu, Int. J. Mod. Phys. A19, 3439 (2004) and its references. 

5. T. Nakano and K. Hicks, Mod. Phys. Lett. A19, 645 (2004). 

6. A.W. Thomas, K. Hicks and A. Hosaka, Prog. Theor. Phys. Ill, 291 (2004). 

7. F. Csikor, Z. Fodor, S.D. Katz and T.G. Kovacs, J HEP 0311, 070 (2003). 

8. S. Sasaki, Phys. Rev. Lett. 93, 152001 (2001). 

9. T.W. Chiu and T.H. Hsieh, hep-ph/0403020. 

10. N.Mathur, F.X.Lee, A.Alexandru, C.Bennhold, Y.Chen, S.J.Dong, 
T.Draper, I.Horvath, K.F.Liu, S.Tamhankar and J.B.Zang, hep-ph/0406196. 

11. H. Suganuma, T.T. Takahashi, F. Okiharu and H. Ichie, Proc. of QCD Down 
Under, Adelaide, March 2004, Nucl. Phys. B (Proc. Suppl.) in press; F. Ok- 
iharu, H. Suganuma and T.T. Takahashi, hep-lat/0407001. 

12. T.R. Klassen, Nucl. Phys. B533, 557 (1998). 

13. H. Matsufuru, T. Onogi and T. Umeda, Phys. Rev. D64, 114503 (2001). 

14. Y. Nemoto, N. Nakajima, H. Matsufuru and H. Suganuma, Phys. Rev. D68, 
094505 (2003). 

15. N. Ishii, H. Suganuma and H. Matsufuru, Phys. Rev. D66, 094506 (2002); 
Phys. Rev. D66, 014507 (2002). 

16. N. Ishii, T. Doi, H. Iida, M. Oka, F. Okiharu, H. Suganuma, hep-lat/0408030. 

17. I. Montvay and G. Minister, "Quantum Fields on a Lattice", (Cambridge 
Univ. Press, Cambridge, England, 1994), p. 1. 



